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1 Resumo 

Neste documento é apresentado o plano de projeto para estágio pós-doutoral a 
ser vinculado ao projeto Pró-Defesa (CAPES - Coordena§ao de Aperfei§oamento 
de Pessoal de Nivel Superior). As atividades partem do estudo de compressao de 
dados SAR brutos antes de processamento radar para síntese de imagens, como 
urna alternativa á compressao da imagem sintetizada após o processamento radar. 
Neste processo, a compressao garante o melhor desempenho para transmissáo e 
armazenamento de dados massivos em opera 96 es embarcadas, além de viabilizar 
a capacidade de opera 9 ao em tempo-real no processamento radar. Na proposta 
aqui descrita, a caracteriza 9 ao de alvo no dominio do sinal bruto é explorada no 
método de análise de transformada. Dessa forma, é possivel alcan 9 ar ganhos de 
eficiencia de compressao e transmissáo e a redu 9 áo de complexidade no proces¬ 
samento radar, dada a garantia de fidelidade ñas regioes de interesse do alvo. 
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2 Introdu^áo e Justificativa 

o radar de abertura sintética {Synthetic aperture radar - SAR) é um sistema de 
imageamento ativo bascado no processamento coerente de múltiplos sinais de eco 
capturados no percurso de urna plataforma em movimento (aeronave ou satélite), 
sendo amplamente utilizado em aplica§ 5 es civis e militares. O SAR transmite 
sinais chirp e os ecos recebidos sao amostrados dentro de componentes In-phase 
(I) e Quadratura (Q), referenciados como dados SAR brutos JIl-O- Sistemas 
SAR modernos produzem imagens de resolu96es muito altas a partir de um grande 
volume de dados brutos processado^ fazendo com que a compressao do sinal 
bruto seja um requisito essencial para redu9ao do volume de dados armazenados 
onboard, a transmissao eficiente para a esta9ao de solo e capacidade de opera9áo 
em tempo-real. Isto pode ser alcan9ado com técnicas de compressao aplicadas ao 
sinal SAR bruto, que normalmente exploram sua característica de ruido devido 
ao espalhamento adicionado incoerentemente com fase e amplitude desconheci- 
das a-Q. 

As duas principáis abordagens de compressao de dados SAR empregam téc¬ 
nicas no dominio espacial e no dominio de transformada. No dominio espacial, 
sao aplicados métodos de quantiza9ao diretamente ao sinal, tais como o BAQ 
(Black Adaptive Quantization) O, BFPQ (blockfloating-pointquantizer), BAVQ 
(block adaptive vector quantization), VQ (vector Quantization) e TCQ (Trellis 
Coded Quantization), sendo o BAQ mais amplamente reconhecido método de 
compressao de dados SAR brutos, bascado na quantiza9áo escalar controlada 
pelas estatísticas de blocos do sinal SAR. No entanto, para garantir melhor desem- 
penho de compressao e qualidade, as técnicas de codifica9ao no dominio da trans¬ 
formada sao consideradas, onde as dependencias estatísticas existentes no sinal 
original podem ser melhor exploradas por urna transforma9áo otimizada seguida 
de quantiza9áo. Dentro dessa abordagem, as técnicas aplicadas a sinais SAR bru¬ 
tos sao bascadas na FFT (fast Fourier transform), DCT (discrete cosine transform) 
e wavelets [O . 

Neste trabalho, a opera9áo radar aplicada a detec9ao de alvos é considerada, 
admitindo um grau de degrada9ao ñas regiSes complementares ao alvo. Dessa 
forma, um compromisso entre a alta resolu9ao na regiao de alvo, através da avali- 
a9áo da imagem sintetizada, e o mais alto grau de redu9ao do dado comprimido 
sao considerados para o desenvolvimento do método de compressao. Dentro desse 
contexto, métodos bascados em transformadas permitem a adequada caracteri- 
za9ao de alvos no dominio SAR bruto, consequentemente, a redu9ao do range 

'Para produzir imagens coerentes a partir de dados SAR brutos, algoritmos de processamento 
radar sao aplicados. No dominio da frequéncia, algoritmos para imageamento tais como Range 
Doppler, chirp scaling e wfc sao amplamente empregados. No dominio temporal, o algoritmo 
backprojection m oferece eficiencia para opera^áo em tempo-real associada a mais alta fidelidade. 
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dinámico médio do sinal SAR resulta no menor tempo de processamento radar, 
eomo um requisito para cenários para opera§áo em tempo-real, o que pode ser 
avallado sob diferentes proeessadores aplieados. 


3 Objetivos 

o objetivo desta pesquisa de pós-doutorado é propor urna caracteriza§áo da 
regiáo de alvo diretamente no sinal SAR bruto no dominio da transformada e, 
adequadamente quantizar para preservar a melhor resolu9áo nesta regiáo de inte- 
resse, eonforme estatístieas do sinal. Consequentemente, a análise de desempenho 
de imagens sintetizadas e a eficiéneia da opera9áo do processamento radar, em 
termos de tempo e eomplexidade, devem ser tais que o opera9áo de imageamento 
e gera9áo das imagens deve oeorrer em tempo-real eom um período mínimo de 
laténeia que será definido pelos requisitos de opera9áo do sistema SAR. Dentro 
deste contexto, busea-se a melhor performance de transmissáo em um cenário de 
sensoriamento remoto e a melhor qualidade das imagens reeuperadas para opera- 
95es de pós-processamento empregadas dentro do interesse do projeto. 

Dessa forma, os objetivos espeeífieos considerados no plano do pós-doutarado 
sao: 

1 . Caraeteriza9áo de alvos no dominio da transformada de sinais SAR brutos; 

2 . Compressáo baseada em quantiza9áo ponderada as regioes do alvo carae- 
terizado; e 

3 . Adapta9áo do sinal reeuperado ao algoritmo de proeessamento radar sele- 
eionado. 


4 Pesquisas Iniciáis 

Considerando os objetivos apresentados na Se9áo[^ as pesquisas inicialmente 
realizadas destaeam a fundamenta9áo de compressáo de dados SAR brutos através 
das prineipais téenicas de compressáo apresentadas na literatura, dentre as quais, 
o BAQ e W-BAQ. Em particular, o método bascado em análise multiresolu9áo do 
W-BAQ é destacado para a caracteriza9áo do alvo no dominio da transformada 
proposto neste plano. A partir da revisáo apresentada ñas próximas se9oes, sáo 
desenvolvidas as metodologías sugeridas na proposta de projeto da Se9áo[^ 
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4.1 Compressáo de Dados SAR Brutos no Dominio do Tempo 

Em muitas aplica§6es de compressáo com perdas existe a necessidade de re¬ 
presentar as saídas da fonte usando um número pequeño de palavras-código, urna 
vez que o número de valores distintos da saída da fonte sao muito maiores do que 
as palavras-códigos disponíveis para representá-las. A quantiza9ao permite repre¬ 
sentar esses infinitos valores em um conjunto de dados muito menores. Usando 
métodos de quantiza9ao adaptativa é possível alcan9ar urna efetiva compressáo de 
dados para urna fonte com parámetros variáveis no tempo. 

Nesse método, as estatísticas da fonte sao estimadas e ajustas para a opera9ao 
otimizada do quantizador. Como a quantiza9áo causa muitas perdas de infor- 
ma9áo, o valor quantizado reconstruido é diferente do valor do sinal de entrada 
do quantizador. Tal diferen9a é chamada de erro de quantiza9áo, medido através 
da teoria de taxa-distor9áo por d{x, Q{x)), representando a distor9áo instantánea 
de Q na entrada x. Para urna variável aleatoria X, a distor9ao esperada é dada por: 


D{Q) = Ed{X,Q{X)). ( 1 ) 

Nesse caso, a fonte X tem fun9ao densidade de probabilidade (PDF) f{x) e, pode 
ser reescrita como: 


Ed{X,Q{X)) = = / d{x,Q{x))f{x)dx 


N 


^ / d{x,Q{x))f{x)dx 

i=l 


rxi 

= / d{x,Q{x))f{x)dx. 

i=i 


( 2 ) 


Assumindo-se que d{x, y) = {x—y)‘^, entáo D{Q) é chamado erro médio quadrático 
(MSE) át Q Q, e = Q{X) — X é entáo definido como o erro de quantiza9áo, desde 
que 


Q{X) = X + Q{X) -X = X + e ( 3 ) 

D{Q{X)) = E[{X - Q{X)f] = E{e^) ( 4 ) 

Outra medida de qualidade derivada da distor9áo de quantiza9áo é a rela9áo 
sinal-ruído (SNR), representada por: 
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( 5 ) 


SNR = 101og.„§^ dB. 
ou comparando com a potencia do sinal de entrada, 

SNRdB = 10 log , (6) 

onde é a variáneia do sinal da fonte. 

É possível observar das rela§6es apresentadas que, menores intervalos de de- 
eisao terao menores erros (£ — e, eonsequentemente, MSE de quantiza§ao a^. 
Dessa forma, O posieionamento ótimo dos limites de deeisao e valores quanti- 
zados sao atingidos com minimiza9áo de cr^ em quantizadores escalares étimos, 
definidos eomo: 


Ed{X,Q*{X)) = min Ed{X,Q{X)), ( 7 ) 

Q&Qn 

onde Qn denota de todos iV-níveis quantizadores e Q* G Qn é o quantizador 
ótimo. 

4.1.1 O Algoritmo Block Adaptive Quantization - BAQ 

O diagrama de blocos do algoritmo BAQ é apresentado na Figura Neste 
esquema, o quantizador Lloyd-Max Q representa o principal módulo do sistema 
e possui urna adapta9ao para as estatístieas bloco-a-bloco do sinal. 

Básicamente, o algoritmo consiste em (i) representar o dado SAR eomplexo 
em seus eomponentes eartesianos (I e Q), (ii) dividir cada cada canal em bloeos 
relativamente pequeños e (iii) adapta9áo dos parámetros de quantiza9áo dentro de 
eada bloeo (threshold e valores de reeonstru9áo) para utiliza9áo no quantizador 
Lloyd-Max, dessa forma, reduzindo o número requerido de bits para eada bloeo. 
Os detalhes do algoritmo sao apresentados nos passos a seguir. 


1 . Estima9áo da variáneia do bloeo; 

2 . Cálenlo do threshold ótimo; 

3 . Quantiza9áo do dado usando os thresholds ótimos eom 2 , 3 , 4 , 5 , 6 ou 8 bits; 
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BAQ 



Figura 1 : Esquema de eodifiea§ao dos dados SAR brutos realizada pelo algoritmo 
BAQ a. 

4 . Transmissao do dado quantizado eom a varáneia; 

5 . Cáleulo dos níveis de trheshold étimos na decodifioa§áo; 

6 . Reoonstru§ao do dado quantizado. 

O quantizador é controlado pelas estatísticas do sinal de entrada que podem ser 
especificadas por um único parámetro de sua distribuÍ9áo de probabilidade (por 
assumir ser gaussiana). Esse parámetro é selecionado para caracterizar a magni- 
tude média do sinal, que pode descrita por 


e 


iV-l 

|/| = |Q| = 2 5^(X„ + -) / 

n=0 




p{x)dx 


( 8 ) 


p{x) 


1 

- p 2cr2 

y/2ña ’ 


( 9 ) 


onde a; é a máxima magnitude média do sinal e a é o desvio padráo do sinal. Po- 
dendo ser simplificado para: 


I\ 


\ Q \ = 127.5 


127 

E- 

n=0 


'(S) 


( 10 ) 
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Estruturalmente, no esquema do BAQ, as amostras do sinal SAR sao divididas 
em blocos para estima9ao das estatísticas e quantiza9ao que, é otimizada para cada 
bloco. Dessa forma, o desvio padrao de cada bloco é estimado, a partir de 


a = 


\ 


MN 


M-l N-l 

EZ 

i=0 j=0 


^2 


( 11 ) 


onde M, N sao os números de amostras das linhas e colunas de cada bloco, re¬ 
spectivamente. Para urna opera9ao ótima, os valores de threshold sao propor¬ 
cionáis ao desvio padrao de cada bloco, sendo esses valores determinados pelo 
cálculo da magnitude média de um bloco usando valores tabelados. Dessa forma, 
a quantiza9áo é feita pela compara9áo de urna amostra da entrada com o conjunto 
de thresholds. Assim, para cada bloco, os níveis de saída do quantizador nao- 
uniforme sao multiplicados pelo desvio padrao do bloco e, na descompressáo, o 
mesmo desvio padrao de cada bloco é transmitido e usado para a reconstru9áo do 
sinal. 


4.2 Compressáo de Dados SAR Brutos no Dominio da Trans¬ 
formada 

A transformada é um método que decompoe um sinal em termos de coefi¬ 
cientes descorrelacionados que podem ser quantizados independentemente de ou- 
tros, demandando diferentes aloca96es de bits de quantiza9ao conforme o grau de 
importancia. A melhor base de sinal para a decomposÍ9áo, i.e., a com melhor com- 
pacta9áo da energía do sinal dentro de poneos coeficientes apresenta a mais alta 
taxa de compressáo para urna dada fidelidade. Quando os coeficientes sao gaus- 
sianos, a base de Karhunen-Loéve (KLT) é ótima, no entanto, a alta complexidade 
computacional da KLT é inviável para aplica95es sobre dados SAR brutos. Dessa 
forma, outras bases ortonormais para decomposÍ9áo do sinal sao aplicadas, tais 
como Transformada de Fourier (TF), a transformada discreta do cosseno (DCT) e 
a transformda wavelet (DWT). 

Específicamente, a DWT é bascada no principio de que a descorrela9áo efi¬ 
ciente pode ser alcan9ada por splitting o dado dentro de duas subsequéncias con- 
tendo informa95es de passa-baixas e passa-altas do sinal original, sendo que a 
maior energía do sinal é carregada pela sequéncia passa-baixas, que garante a 
compacta9áo de energía necessária. Na aplica9áo para compressáo de dados SAR 
o uso de transformadas sáo bascadas na normaliza9áo de blocos dos coeficientes 
I e Q, que assumem o modelo de um processo Gaussiano estacionário com den- 
sidade espectral de potencia (psd) ñas dire95es range e azimute. Particularmente, 
a psd na dire9áo range é definida pelo comprimento de banda do pulso chirp, en- 
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quanto que a psd na dire§áo azimute é proporcional ao padrao da antena azimute. 
Essas pdf’s apresentam leve comportamento passa-baixas, mais evidente na di- 
re§ao range, demonstrando a existencia de certa quantidade de correla§ao dentre 
as amostras que podem ser exploradas pelo principio da transformada [Hl . 

O método apresentado na Figura [^descreve o processo de compressáo de da¬ 
dos SAR com a aplica§áo de transformada. Considerando que os coeficientes da 
transformada sao distribuidos gaussianamente, o BAQ pode ser aplicado para a 
quantiza9áo dos dados. Dessa forma, os coeficientes da transformada de cada 
bloco sao quantizados e transmitidos com o desvio padrao calculado. 


BAQ 



Figura 2 : Esquema de codifica9ao dos dados SAR brutos no domñio da trasfor¬ 
mada (DCT/wave/et) seguida da quantiza9ao realizada pelo algoritmo BAQ. 


4.3 Medidas de Qualidade 

Os parámetros de básicos para avalia9áo de qualidade e compara9áo entre al¬ 
goritmos de compressáo de sinais SAR complexos sao a taxa de compressáo (CR), 
rela9áo sinal-ruido (SNR), Erro Médio de Fase (MPE) e o Desvio Padráo da Fase 
(PSD), além de medidas de qualidade observadas ñas imagens sintetizadas gera- 
das a partir do sinal SAR descomprimido. As formula96es desses parámetros sáo 
descritas a seguir: 


SNR 


lOlog 


1 sr^Ka sr^Kr 2 
KaKr 2^i=\ l^j=l -^ij 


KaKr = 1 


.X. 


V 


X. 


tj) 


( 12 ) 
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MPE 


1 


(13) 


KaKr 


Ka Kr 

\^ij ~ 

i=l j=l 


Onde Xij e Xij sao as amostras do dado e o dado codificado, respectivamente. 
Ka e Kr sao o número elemento em azimute e range. 


PSD 


1 



Ka Kr 


i = l j = l 




( 14 ) 


A MPE e o (PSD) sao métricas importantes para aplica§ 5 es interferométricas 
em que a acuracidade das informa96es de fase do sinal é essencial para a gera§áo 
adequada das imagens. 

Um exemplo de aplica9ao com avalia9áo da compressao do sinal SAR bruto 
submetido ao algoritmo BAQ é apresentado na Eigura e na Tabela Neste 
exemplo, a simula9áo de desempenho do BAQ é realizada usando software IDE e 
Matlab e, empregando dados SAR complexos brutos, onde os canais I eQ do sinal 
sao submetidos individualmente ao algoritmo BAQ em blocos de 32 x 32 para 
extra9ao das medidas estatísticas aplicadas ao quantizador. Como resultado, o 
sinal SAR recuperado e submetido ao processamento radar é capaz de reconstruir 
a imagem com baixa degrada9ao, conforme observado na Eigura 


Tabela 1 : Desempenho do Algoritmo BAQ usando 4 bps aplicados aos canais 
SAR no dominio do sinal bruto em amostras de 8192 x 2048 em range e azimute 
cedidas pela BRADAR. 


Canal SAR 

CR 

SNR 

MPE 

PSD 

/ 

2.7 

14.0 

11.56 

10.0 


2.7 

13.9 

11.34 

10.2 
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(a) (b) 


Figura 3: Imagens geradas pelo processamento SAR Range Doppler a partir do 
sinal SAR complexo (dados brutos) em cena de 8192 x 2048 amostras em range 
e azimute cedidas pela BRADAR: (a) Imagem sintetizada a partir do sinal SAR 
original e (b) Imagem sintetizada a partir do sinal SAR recuperado do BAQ com 
blocos de 32 x 32. 

5 Proposta de Projeto 

A proposta de pesquisa aborda a transmissáo eficiente de sinais SAR bru¬ 
tos, obtidas de dataset específico do projeto, através de técnicas de compressao e, 
ainda, a redu 9 ao observada no tempo de processamento da síntese de imagens no 
processamento radar com os dados SAR brutos recuperados. Particularmente, sao 
considerados para a compressao a abordagem BAQ no dominio da transformada 
wavelet. Dentro deste contexto, busca-se desenvolver um método de carateriza§ao 
das regioes de alvos potenciáis para aplica§ao em urna análise multiresolu§ao no 
dominio do sinal SAR bruto. Urna abordagem previamente realizada pelo autor 
em dH aplica técnica de denoising em sianis SAR brutos, porém, sem a caracte- 
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riza§ao de alvos. A correta descri§áo do alvo seleciona os coeficientes da regiáo 
de interesse para urna quantiza9ao de ordem superior, mantendo menores níveis 
de quantiza§ao ñas regioes complementares. Neste ponto, prop6e-se um limiar 
de sele§áo de níveis de quantiza§ao ponderados aos coeficientes correspondentes 
as regioes de alvos potenciáis. Dessa forma, a regiao mais provável do alvo de 
interesse é mantida com a mais alta resolu§ao e melhores índices de quantiza§ao 
no processo de codifica§ao, o que garante maiores taxas de compressao. 

Urna consequéncia direta é a simplifica9áo no processamento radar, com um 
menor volume de dados a serem processados, dessa forma, considerando os pro- 
cessadores range-doppler e Backprojection como referencias, pretende-se obser¬ 
var a capacidade para opera9ao de tempo reduzida em rela9ao ao sinal SAR bruto 
nao comprimido. Tal análise parte das medidas de tempo dentro dos intervalos de 
tempo de compressao e processamento radar, descritos na Figura]^ a seguir. 



Com Compressao 

AT' 


Sem Compressao 

AT 


Figura 4 : Análise de opera9áo em tempo real do sistema SAR considerando o 
fluxo do sinal SAR bruto para síntese e os tempos de codifica9áo e processamento 
envolvidos. 

Consequentemente, espera-se obter: 

1 . Caracteriza9áo potencial do alvo no dominio da transformada; 

2 . Método de quantiza9áo ponderada aos coeficientes do sinal SAR transfor¬ 
mado potenciáis de alvo (regiáo do alvo) 

3 . Aplica9áo no processamento radar para opera9áo em tempo-real (proces- 
sador temporal e range doppler) 

4 . Fidelidade da imagem de síntese para aplica96es de interesse (detec9áo de 
alvo, classifica9ao, etc.) 
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6 Cronograma de Atividades Previstas 

De urna forma geral, o cronograma apresentado na Tabela|^descreve as ativi¬ 
dades previstas para o estágio pós-doutoral, compreendendo, inicialmente, um 
período de 7 meses (maio/2019 a dezembro/2019). 


Período 

Atividades 

Mai. 2019 

• Plano de Trabalho do Pós-Doutorado e Documenta 9 áo para implementa 9 áo 
da bolsa. 

• Integra 9 áo com o grupo de pesquisa e pesquisas iniciáis. 

Jun. 2019 

• Aplica 9 ao de métodos de compressao de sinais SAR brutos em conjuntos de 
dados disponibilizados para o Projeto. 

Jul. 2019 

• Caracteriza 9 ao do alvo no sinal SAR bruto no dominio da transformada —> 
Resultado # 1. 

Ago. 2019 

• Método de decisao para quantiza 9 ao em fun 9 áo da regiáo do alvo caracteri¬ 
zada no dominio da transformada (sinal SAR bruto) —)> Resultado # 2. 

Set. 2019 

• Análise de desempenho do sinal comprimido e recuperado após o processa- 
mento radar (imagem SLC). 

Out. 2019 

• Simula 96 es e testes para o desempenho de opera 9 ao em tempo-real —í- tempo 
de codifica 9 áo/decodifica 9 áo -|- tempo de processamento radar. 

Nov. 2019 

• Elabora 9 ao de abstract estendido envolvendo o resultado # 1 para confer¬ 
encia internacional IGARSS (International Geoscience and Remóte Sensing 
Symposium)/2020. 

• Elabora 9 ao de relatório do pós-doutorado. 

Dez. 2019 

• Elabora 9 ao de artigo envolvendo o resultado # 2 em periódico das sociedades 
IEEE Geoscience and Remóte Sensing (GRSS), IEEE Signal Processing Soci- 
ety (SPS) ou IEEE Communications Society (ComSoc). 

• Eechamento do relatório do pós-doutorado. 


Tabela 2: Cronograma de atividades previstas para o estágio de pós-doutorado 
considerando um período inicial de 7 meses. 
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7 Considera^oes Fináis 

o método de compressáo de sinais SAR brutos orientada ao alvo proposto 
neste plano considera a necessidade de transmissao eficiente e opera 9 ao em tempo- 
real de imagens SAR sintetizadas, o que viabiliza a implementa 9 áo de proces- 
sadores radar com tempo reduzido e com qualidade das imagens de síntese ad- 
equadas por manter alta resolu 9 áo ñas regioes de interesse, reduzindo a com- 
plexidade ñas análises de pós-processamento, tais como a detec 9 ao de alvos e 
classifica 9 ao. Adicionalmente, os desdobramentos da pesquisa proposta incluem 
a integra 9 ao com pesquisas relacionadas dentro do grupo de trabalho do pro¬ 
jeto Pró-Defesa e a estrutura 9 ao de laboratorio para simula 9 ao e ensaios para 
prosseguimento de futuras pesquisas na área de radares no contexto académico- 
institucional. 
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